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We study the inertial modes of slowly rotating, fully relativistic compact stars. The equations 
04 ' that govern perturbations of both barotropic and non-barotropic models are discussed, but we 

, present numerical results only for the barotropic case. For barotropic stars all inertial modes are 

' a hybrid mixture of axial and polar perturbations. We use a spectral method to solve for such 

modes of various polytropic models. Our main attention is on modes that can be driven unstable 
by the emission of gravitational waves. Hence, we calculate the gravitational-wave growth timescale 
for these unstable modes and compare the results to previous estimates obtained in Newtonian 
gravity (i.e. using post-Newtonian radiation formulas). We find that the inertial modes are slightly 
stabilized by relativistic effects, but that previous conclusions concerning eg. the unstable r-modes 
remain essentially unaltered when the problem is studied in full general relativity. 
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^ . I. INTRODUCTION 

<N . 

That gravitational waves can drive various modes of oscillation in a rotating neutron star unstable was first suggested 
by Chandrasekhar The detailed mechanism behind this instability was explained by Friedman and Schutz |2[ ^, 
who showed that the instability sets in when an originally retrograde mode becomes prograde (according to an inertial 
' observer) due to the rotation of the star. They also showed that the instability is generic, that, for arbitrarily slow 
r. rotation, any perfect-fluid stellar model has unstable modes with sufficiently large values of the azimuthal eigenvalue 
' m (the mode depends on if as exp(im(/3)). The existence of this radiation-driven instability is potentially important 
qh, since it could limit the attainable spin rate of astrophysical neutron stars Q|. However, a detailed assessment of its 
^ ' astrophysical relevance is complicated by the fact that viscosity tends to counteract the growth of an unstable mode. 
One must account not only for the familiar hydrodynamic shear and bulk viscosities 1 1] , but also exotic mechanisms like 
the mutual friction that is relevant in superfluid neutron stars. Once a star has cooled below the superfluid transition 
temperature, mutual friction appears to suppress the instability of f- modes , and gravitational waves from f- modes 
appeared to set a limit on rotation only slightly more stringent that the maximum spin of an equilibrium model. 
^ ' Recent work has modified this conclusion in two ways. First, numerical studies of the marginally stable "neutral" 
f-modes of fully relativistic, rapidly rotating polytropes by Stergioulas and Friedman Q showed that relativistic effects 
destabilize the modes of a rotating star considerably. Most importantly, one finds that in general relativity the m = 2 
f-mode may have a neutral point for attainable rates of rotation. This is in contrast with the Newtonian result that 
the m = 2 mode is unlikely to become unstable in uniformly rotating stars and it is important since the quadrupole 
mode is the most efficient emitter of gravitational radiation. Hence one would expect it to lead to the fastest growing 
instability. For realistic equations of state it has been shown ||] that in a typical I.AMq star the m — 2 f-mode has a 
neutral point near « 0.85Qk where ^Ik represents the mass-shedding limit. 

The second piece of evidence follows from the fact that the inertial r-modes are unstable at any rate of rotation in a 
perfect fluid star The great surprise of a few years ago was the discovery that, despite the fact that they radiate 

mainly through the current multipoles, the unstable r-modes could potentially limit the spin of a rotating neutron 
star significantly Since its original discovery the r-mode instability has been discussed in a number of papers, 

and we refer the interested reader to recent review articles for detailed discussions of the literature [ O , [l 4 | p^ , p^ , p7[ . 
At the present time key issues concern the nonlinear evolution of an unstable mode |18, |l^ |2]], p2[, the possible 



presence of hyperons in the neutron star core (since hypcrons may lead to a very strong bulk viscosity) |2^, |2J^ 
and the role of superfluidity ^ . 

The present paper concerns the effect that general relativity has on the instability of the r-modes and other inertial 
modes. Intuitively, one would expect this to be a relevant issue since the relativistic framedragging is an order 
effect which would affect inertial modes at leading order. In the last couple of years significant progress toward an 
understanding of the nature of inertial modes in general relativity has been made. The picture that is emerging is 
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largely similar to that in Newtonian gravity. In the Newtonian case the r-modes are a purely axial parity subset 
of the large class of inertial modes which generally have velocity fields that are described by a hybrid mixture of 
axial and polar components to leading order in fi. Non-barotropic Newtonian stars have an infinite set of r-modes 
for each permissible I and m (corresponding to the spherical harmonic Y)™ used to describe the velocity field), while 
barotropic models retain only a vestigial r-mode for each I — m, see |2^, for detailed discussions. In relativity one 
can prove that barotropic stars (where the perturbations are described by the same one-parameter equation of state 
as the background star [p9|) have no purely axial inertial modes pO[ |. All inertial modes of such stars are hybrids. 

As in Newtonian theory the relativistic non-barotropic r-mode problem is significantly different. Like their 
Newtonian counterparts, relativistic non-barotropic stars have modes that are purely axial in the spherical limit. These 
modes are determined by solving a single ordinary differential equation for one of the perturbed metric components , 
cf. eq (^). We have previously shown |^ that discrete mode-solutions to this equation exist for uniform density 
stars. For more realistic equations of state, the problem is complicated by the fact that this equation corresponds to a 
singular eigenvalue problem which also admits a continuous spectrum ||3l|, ^ . The dynamical role of this continuous 
spectrum, or indeed if it remains present when more physics is included in the model, is not clear at the present 
time. In order to determine a purely axial mode of a non-barotropic star one must typically identify a discrete mode 
embedded in the continuous spectrum. This technical difficulty has led to suggestions that the r-modes may not even 
"exist" for certain relativistic stars jH, 0. However, from eq. (|2^), it is clear that the slow-rotation approximation 
is no longer consistent in regions where a — w ~ 0{fl^) or smaller. This means that the problem likely requires a 
"boundary layer" approach ||35| where either f2^ terms, viscosity or the coupling to polar perturbations are included 
in the analysis of the region near the singular point. Recent results by Yoshida and Lee and Ruoff et al. 
support this view. 

So far there have been two studies of the growth timescale for the unstable relativistic r-modes |3^, Both 
concern the non-barotropic problem. The results suggest that post-Newtonian estimates of the instability growth 
time are surprisingly good, but also indicate a weakening of the instability once the star becomes very compact. 
The aim of the present paper is to extend these results by considering the effect of radiation reaction on various 
inertial modes of barotropic stellar models. We expect the results obtained from a barotropic model to be relevant 
even though real neutron stars will have internal stratification associated with composition gradients (eg. due to the 
varying proton fraction) . As discussed by Reisenegger and Goldreich |^ this will lead to a Brunt- Vaisala frequency 
N ~ 500 s~^{p/ Pnucicar)^^^ and the presence of core g-modes with frequencies in the range 100-200 Hz. In a very slowly 
rotating star the buoyancy force will dominate the Coriolis force and the g-modes will remain largely unchanged, but 
when >> N the situation will be reversed and the g-modes will be almost entirely rotationally restored. In other 
words, one would expect all low frequency modes to be well described by the inertial modes of a barotropic model 
for rotation periods shorter than 2-3 ms. This means that the results obtained in this paper should be a reasonable 
approximation for neutron stars in Low-Mass X-ray binaries (expected to have rotation rates in the range 250-500 Hz), 
and an accurate representation of the inertial modes of a newly born neutron star spinning at or near the break-up 
limit. 

The layout of the paper is as follows: In section II we summarize the equations that describe the oscillations of 
a slowly rotating relativistic star. Section HI provides a description of the numerical method we use to solve the 
eigenvalue problem, while the obtained results are discussed in Section IV. Section V is devoted to a discussion of our 
method for estimating the radiation reaction timescale for inertial modes. Our conclusions are given in Section VI. 

Throughout the paper we use the following conventions: We refer to ref. as Paper I, ref. as Paper II , while 
Paper HI is | ]35[ |. Our numbering convention for equations from these papers is such that e.g. (11,2.5) means Eqn. 
(2.5) from Paper II. Unless otherwise stated we use geometrized units c = G = 1. 

II. PERTURBATIONS OF SLOWLY ROTATING STARS 
A. The Equilibrium Star 

We consider a perfect fluid star rotating slowly with uniform angular velocity Q. By slow rotation we mean the 
assumption that Q is small compared to the Kepler velocity, fix — 0.67-\/7rGe, at which the star is unstable to mass 
shedding at its equator (e represents the average energy density in the star). In particular, we neglect all quantities 
of order or higher. In this approximation the star retains its spherical shape, because the centrifugal deformation 
of its figure is an order Q'^ effect . The only new order effect that arises because of general relativity is the 
rotational framedragging, denoted by i-u{r) below. 

To first order in Q the equilibrium state is described by a stationary, axisymmetric spacetime with metric, gap, of 
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the form ^ || 

ds^ = -e^^^^Ut^ + e^^^^'Ur^ + r^dO"^ + r^sin^Odip^ - 2uj{r)r^ sin^Odtdip (1) 
with perfect fluid matter source 

Tap = (e + p)uaup + pgap- (2) 

Here, e and p are, respectively, the total energy density and pressure of the fluid as measured by an observer moving 
with unit 4-velocity, 

= e-''(i" + J7</j"); (3) 

— {dt)" and = being, respectively, the timelike and rotational Killing vectors of the spacetime. 

The metric and fluid variables are required to satisfy Einstein's equation, GafS = SirTap. This reduces to the 
well-known Tolman-Oppenhcimer-Volkov (TOV) equations for the "spherical quantities" (those of order fi^) together 
with Hartle's equation for the order quantity, 

ij{r) = Lj{r), (4) 

that governs the dragging of inertial frames induced by the rotation of the star. For the exact form of these equations, 
we refer the reader to Paper II; Eqs. (11,3.4-3.7) and (11,4.3). 

To complete our specification of the equilibrium star we must provide an equation of state (EOS) relating the density 
and pressure. In this paper, we will always require our equilibrium solution to satisfy a one-parameter EOS, e — e{p). 
This is an accurate assumption for equilibrium neutron stars since their temperature is likely to be significantly below 
the Fermi temperature Tp ^ 10^^ K. For simplicity we use the polytropic EOS, 

p = Kp'+i 

e = p + np, (5) 

where p is the rest-mass density, n is the polytropic index and K is the polytropic constant. We use a set of polytropic 
indices (n — 0.0, 0.5, 1.0, 1.5) that span the range of compressibilities of proposed realistic neutron star equations of 
state g, ||. 

Although our equilibrium model obeys a one-parameter EOS, we do not necessarily require the perturbed fluid 
to satisfy the same EOS. For an adiabatic perturbation of an equilibrium star obeying a one-parameter EOS, the 
perturbed pressure and energy density are customarily related by 

r^-T^,+CAa (6) 

Tip [e+p) 

where ri(r) is the adiabatic index, is the Lagrangian fluid displacement and where the Schwarzschild discriminant, 

Ac = j^-^Vae - -^-V„p, (7) 

[e+p) Tip 

governs convective stability in the star. In general, the adiabatic index Fi need not be equal to the constant 

r^ii±^^^i-Hi (8) 

p de n 

associated with the equilibrium (polytropic) EOS. In terms of this constant we have 

A. ^ (1 - i-) iv.p. (9) 

We will call a model barotropic if and only if the perturbed configuration satisfies the same one-parameter EOS as the 
unperturbed configuration. In this case Fi = F and the Schwarzschild discriminant vanishes identically. Such stars 
are marginally stable to convection, and since they have no internal stratification they do not admit finite frequency 
modes restored by buoyancy (g-modcs) pO[ . In this paper we will consider perturbations of both barotropic and 
nonbarotropic stars, but will focus mainly on the barotropic case in our numerical work. The non-barotropic case has 
already been discussed in detail in [pOl Esl, p5|. 
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B. The Perturbation Equations 



We now consider the rotationally restored (inertial) modes of a slowly rotating relativistic star. The equations 
governing such perturbations were derived in detail in Paper II. Here we will simply quote the main results needed 
for the present analysis. 

Since the equilibrium spacetime is stationary and axisymmetric, we may decompose our perturbations into oscil- 
lation modes proportional to e*('^*+™'^). For convenience, we will always choose m > 0, since the complex conjugate 
of an m < mode with real frequency cr is an m > mode with frequency —a. Note that a is the mode frequency 
measured by an inertial observer at infinity. 

In the Lagrangian perturbation formalism p5| the basic variables are the metric perturbation haji and the 
Lagrangian displacement We begin by expanding these variables in vector and tensor spherical harmonics. The 
Lagrangian displacement vector can be written 



:M^/(r)y,"V" 



(10) 



where we have defined, 



and introduced the "comoving" frequency, 



(11) 



(12) 



The exact form of expression has been chosen for convenience. In particular, we have used the gauge freedom 
inherent in the Lagrangian formalism |47, 48j to set = 0. Note also the chosen relative phase between the terms in 
( p^ ) with polar parity (those with coefhcicnts Wi and Vi) and the terms with axial parity (those with coefficients Ui). 
Working in the Regge- Wheeler gauge p9|, we express our metric perturbation as 



hill/ — 6 



l—7n 



HiA{r)Y- 

symm 
symm. 



2A 



symm 
symm 



r'^Ki{r)Y{^ 

r^sin^eKi{r)Y/^ 



(13) 



Again, note the choice of phase between the polar-parity components (those with coefficients Hqj, Hi i, and Ki) 
and the axial-parity components (those with coefficients Hq i and hu). 

In Paper II, we also found it convenient to make use of the Eulerian perturbation formalism, whose basic variables 
are the metric perturbation haf3 and the perturbed density, 6e, pressure, Sp, and fluid 4- velocity, Su°'. Since they 
represent scalar quantities the Eulerian changes in the density and pressure may be written as 



Sei{r)Yre'' 



(14) 



and 



Sp=J2SMr)Yre^' 



(15) 



respectively, while the Eulerian change in the fluid velocity may be expressed in terms of ^" and haf} as defined above, 

6u°' = q'^p£u£.^ + \u°'u'^u''hp^, (16) 

where g"'^ = g"^ + u^u^ . 

We showed in Paper II that the rotationally restored modes of a slowly rotating star have a fundamentally different 
character depending on whether the star is barotropic or nonbarotropic. The difference pertains to the character of 
the modes in the limit as the star's angular velocity, fi, goes to zero. This is yet another reason why it is appropriate to 
consider the problem within the slow rotation approximation. In particular, we proved in Paper II that a relativistic 
barotrope does not admit distinct classes of r-modes or g-modes (modes whose limit as f2 — s- are purely axial or 
purely polar, respectively). Instead, the generic inertial mode of such a star is a hybrid mixture of axial and polar 
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components to lowest order in f2. In contrast, non-barotropic stars allow distinct g-modes already in the non-rotating 
case and purely axial r- modes may exist at lowest order in f2. To reflect the fundamental difference between the two 
cases, we organized our slow-rotation expansion by requiring our perturbation variables to obey the following ordering 
in powers of Q, 

Ui,ho.i^ 0(1) 



0(1) 

o(r!2) 



barotropic stars 
nonbarotropic stars 



Ho^i,H2,i,Ki,hi^i,Sei,Spi,a ^ 0{Q) . 



(17) 



The perturbation equations may then be grouped in powers of fl and solved order by order. To compare our results 
with the Newtonian r-modes and hybrid modes we need only find the leading order mode solutions; that is, we need 
only find the mode frequency to order fl and the eigenfunctions to 0(1). As discussed in Paper II, these turn out to 
be determined by a subset of the perturbation equations up to first order in Q. 

It is relevant to point out that the equations derived by this procedure will be somewhat different in different 
regions of spacetime (see Table ^. In the "near zone," the region in which or << 1, we will be able to ignore second 
time derivatives, whereas we cannot do this in the "wave zone," the region in which ar » 1. Because the inertial 
modes that we are primarily interested in are restored by the Coriolis force, their frequencies scale with the angular 
velocity of the star, ct ~ $7. For slow rotation, this implies that the near zone extends far away from the star into 
the nonrclativistic region (M/r << 1) and that the wave zone will be located entirely within the nonrelativistic 
region This will be important in Sect. IV when we calculate the energy radiated in gravitational waves and 

the timescales on which gravitational radiation reaction drives the unstable modes. For now, we will focus on the 
equations that are relevant in the near zone, which were derived in Paper II. These will allow us to find the eigenvalues 
and eigenfunctions of the modes that we are interested in. In Sect. IVB we will consider the equations more generally 
when we derive an expression for the radiated energy. The various spacetime regions are depicted schematically in 
Table |. 





relativistic zone, 


M 1 

r 




nonrelativistic zone, — 


<< 1 




r = 


neutron star 


r — 


R 




exterior spacetime 




r — > oo 








near zone, 


ar ~ S7r << 1 




wave 


zone, err ~ fir >> 1 





TABLE L The spatial regions relevant to the relativistic inertial-mode problem. 



The ordering ( p7| ) for the barotropic case is slightly more general than for the nonbarotropic case because the 
polar-parity coefficients Wi, Vi and Hi_i are not assumed to be negligible compared to the axial-parity coefficients Ui 
and ft-o,i- We will retain all of these variables in presenting the relevant equations and then specialize to each of the 
two case s. (B ecause hi^i is an order SI variable, we ignore it and drop the "0" subscript on /iq,;, writing it as hi. Only 
in Sect. IV B| will it be necessary to restore this distinction.) 

As presented in Paper II, the relevant 0(1) equations (those that apply to the perturbed spherical star) are: 
Eq. (11,3.20), 



Eq. (11,3.23), 



167r(e + p) o\ 



p-(i'+A) , 
Vi = , , , [{e+p)e''+\Wi\ . 



l{l + l){e+p) 



(18) 



(19) 
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and Eq. (11,3.22), 



h'; - {v' + X')h[ + 



hi = -{v' + X')Ui 



(20) 



where a prime denotes a derivative with respect to r. Throughout this work wc have used Eq. ( |l8| ) [i.e., Eq. (11,3.20)] 
to ehminate the metric variable Hi^i{r) in favor of Wi{r). 

To close the system of equations we need to retain only two of the equations that arise at 0(il). As discussed in 
Paper II, the relevant pair are the two independent components of Eq. (11,4.16), which enforces the conservation of 
vorticity in constant entropy surfaces [see Eq. (]2^) below]. The [Qip] component of Eq. (11,4.16) leads to Eq. (11,4.53), 



0= [l{l + l)Kn{hi + Ui) -2mQUi] 

+ {l + l)Qi [^dr (r^oje-^-) Wi_i - 2{l - l)CuVi_, 

'-^dr (r2we-2-) Wi+i + 2{l + 2)LdVi+i 



while the [r9] component of Eq. (11,4.16) gives Eq. (11,4.54) 



= 



(e+p) 



(/ - l)Qi ( iSpi^i + ^Wi^^ -{1 + 2)Qi+i (iSpi+i + ^Wi+i 



+Qi {I - l)Kndr {e-^^Vi-i) - 2mdr {uje-^^Vi-i) 

+ (r^cSe-^-) + (/ - l)«17e-2- (^^^{^ - i) e^^Wi^ 

mnVLdr [e-^^ihi + Ui)] + 2dr {Cje-^^Ui) {{I + l)Qf - lQ^i+^) 
+ ^dr (r2cDe-2-) Ui [m2 + /(/ + 1) [Qf^^ + Qf - l)] 



li+i 



{I + 2)Kndr (e-^'^Vi+i) + 2mdr (we-^'^VJ+i 



The constants Qi were defined in Paper II to be 



{l + m){l - m) 



{21- 1)(2Z + 1) 



1/2 



(21) 



(22) 



(23) 



Notice that in writing Eq. ( p2[ ) we have retained (in the first line) the term containing the Schwarzschild discriminant 
Ar from the right hand side of Eq. (11,4.16), 



(24) 



We have done this in order to emphasize the difference between the barotropic and nonbarotropic cases. To find the 
rotational modes to leading order in 17 we need a complete set of equations involving only our 0(1) variables. If is 
not identically zero, the retained term in (E2) brings about a coupling between these variables and the 0{Vl) variable 
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5pi. In a barotropic star, Ar is identically zero, so this coupling between 0(1) and 0{^l) variables vanishes. In this 
case, the five equations (p^-(p2|) do, indeed, involve only the 0(1) variables hi, Ui, Wi, Vi and Hij. These equations, 
therefore, comprise a complete set and fully determine our normal mode eigenvalue problem. In a nonbarotropic star, 
however, Ar ^ 0, and the coupling between 0(1) and 0{n) variables in Eq. (22) does not vanish. In this case, the 
equations (na)-(22) do not involve only the variables hi, Ui, Wi, Vi and Hi i. There are then two options: The first is 
reminiscent of the situation for r- modes of non-barotropic Newtonian stars. The problem would become well posed if 
we extended the analysis to one order higher in Q and thus determined also Spi etcetera. The second possibility 
would be to obtain a well-defined eigenvalue problem by assuming the non-barotropic ordering ( p^ ) in which only 
the axial variables Ui and hi are 0(1). With this choice, one obtains an eigenvalue problem from Eqs. ( pO| ) and 
( pT| ) - dropping Wi and Vi from the latter as O(ri^) quantities. These equations can be reexpressed as Eq. (11,4.24): 
Kojima's pl| eigenvalue equation governing hi. 



(a 









+ AM 




dr 







hi } + 16TT{e+p)ahi = , 



(25) 



where uj = lu/Q and a = l{l + l)K,/2m and where Ui is determined (once the eigenvalue, a, and hi are known) by Eq. 
(pT|), which becomes. 



Ui = 



(lj — a) 



(26) 



To summarize: The eigenvalue problem governiiig the rotational modes of a barotropic relativistic star to lowest 
order in ft is determined by the set of equations (|l8|)-(|2^) for the variables hi, Ui, Wi, Vi and -ffi,z, which are all 
assumed to be zeroth order in Q. By contrast, the eigenvalue problem governing the rotational modes of a non- 
barotropic star to lowest order in f2 is determined by Eqs. (^5|) and (|26| ) for the variables hi and Ui, which alone are 
assumed to be zeroth order in Q. 

To complete the specification of our eigenvalue problem we must impose appropriate boundary conditions. We 
require our variables to be regular everywhere in the spacetime, which implies that they vanish at the origin, r — 0. 
The fluid variables Ui{r) and Vi{r) are otherwise unconstrained - apart from the fact that they vanish outside the 
star. In Paper II, we derived the boundary condition (11,4.64) on the fluid variable Wi{r) at the surface of the star, 
r = R, 



Wi{R) = {all I). 



(27) 



This condition and Eq. ( [l8| ) imply that the metric variable Hi^i{r) vanishes at the surface, and in the exterior, of the 
star. By inspection, it is clear that all but one of the perturbation equations vanish in the vacuum exterior to the star 
(r > R). The one surviving equation is that governing the metric variable hi[r), i.e., Eq. (|2^) for barotropic stars 
or Eq. ( p5| ) for nonbarotropic stars. In the near zone, where these equations apply, they both reduce to the same 
equation outside the star: Eq. (11,4.66), or. 



(1 



-)■ 



l{l + l) 



AM 



0, 



(28) 



where M is the total gravitational mass of the star. This equation has a regular singular point at r = oo, so it 
has at least one regular series expansion about this point. The other, linearly independent, solution turns out to be 
singular at r = cxd, growing like r'"*"^. As we will see in Sect. IVB wc may ignore this second solution in the near zone 
because it is of order (tTr)^'+^ relative to the nonsingular solution. To flnd the mode eigenvalues and eigenfunctions 
to lowest order in we therefore need only the regular solution. As discussed in Paper II the regular solution can be 
immediately written down as Eq. (11,4.67), or 



Hr)-f:hus if) 



l+s 



where the coefficients hi ^ are given by the recursion relation (11,4.68), 



or, equivalently, by. 



(2M\ (; + s-2)(/ + .s + l)-. 
- \Ji) s{2l + s + l) ^'--^ 



_ {l + s-2)\{l + s + l)\{2l + \)\ / 2M 
~ s\{l ~ 2)1(1 + 1)\{21 + s + 1)1 [it 



H,0 



(29) 



(30) 



(31) 
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with hifi an arbitrary normalization constant. This known exterior solution must be matched to the interior solution 
for hi{r) at the surface of the star. This provides a boundary condition on the interior solution. We require that the 
solutions be continuous at the surface, 

\im[hi{R-e)~hi{R + e)]^0, (32) 

for all I (which fixes the normalization constant /i;,o), and that the Wronskian of the interior and exterior solutions 
vanish at r = i?, i.e. that 

lim [hi{R - s)h'i{R + e)- h'i{R - e)hi{R + e)] = 0, (33) 

e— >0 

for all I. 

Finally, we note that since we are working in linearized perturbation theory there is a scale invariance to the 
equations. If (^"7/10/3) is a solution to the perturbation equations then (5^", S'/Iq,^) is also a solution, for constant 
S. We will sometimes find it convenient to impose the following normalization condition in addition to the boundary 
and matching conditions just discussed: 

UmiT = -R) = 1 for axial-led inertial modes and r- modes, ^.^^^ 
Um+iir = R) = I for polar- led inertial modes. 



III. MODE FREQUENCIES AND EIGENFUNCTIONS 



From the equations in the previous section it is clear that the task of determining the relativistic analogues of the 
Newtonian r-modes is different depending on whether the star is barotropic or non-barotropic. Each problem presents 
it's own computational challenge. For non-barotropic stars one must deal with the fact that Kojima's equation ( [25| ) 
represents a singular eigenvalue problem. The barotropic case is conceptually easier because one does not have to 
worry about singularities. On the other hand, all inertial modes of a barotropic star will have a hybrid nature which 
complicates the numerical solution of the eigenvalue problem considerably. In this paper we focus our attention 
on barotropic stars since: i) one would expect the inertial modes of a more complex stellar model (with internal 
stratification) to be similar to those of a barotropic star for sufficiently rapid rotation, ii) the nonbarotropic case has 
already been discussed by several authors |33, |4, 35 1. 

We have used the spectral method described in Appendix ^ to solve numerically for a sample of rotationally restored 
modes of slowly rotating, fully relativistic barotropes. The relevant set of equations for this problem is Eqs. (pq)-(p2[), 
which comprise a system of ordinary differential equations for the variables Ui{r), V/(r), Wi{r) Hij{r) and hijr) (for 
all allowed I). Together with the boundary and matching conditions at the surface of the star, these equations form 
a nonlinear eigenvalue problem for the parameter k, the dimensionless mode frequency in the rotating frame. [In 
practice, because the metric variable Hi i is related algebraically to the fluid variable Wi by Eq. ([l8|), we eliminate 
Hi I from the other equations and solve the system ([l9|)-(|2^).] 

For simplicity, we have restricted our study to relativistic polytropes even though it would be straightforward to 
generalise our calculation to tabulated realistic equations of state. This is an important step that should eventually 
be taken, but we feel that we should first try to understand the overall effect that general relativity has on the inertial 
modes of a compact star. 

In Newtonian barotropic stars there remained a large set of modes that were purely axial to lowest order in Q: 
the r-modes with spherical harmonic indices I — m. The I ~ m = 2 r-mode is the one expected to dominate the 
gravitational wave-driven instability of sufficiently hot and rapidly rotating neutron stars ||lf|, |l^. We showed in 
Paper II, however, that relativistic barotropes do not admit such modes. Pure r-modes with I = m > 2 are not 
allowed by the perturbation equations (|20|)-(|2^). The corresponding modes are instead axial-led hybrid modes. We 
have previously solved explicitly for these important hybrid modes to first post-Newtonian order in uniform density 
stars [cf. Eqs. (II,5.33)-(II,5.40)]. We now report on a more general numerical study of these and other hybrid inertial 
modes in fully relativistic stars. 

Apart from the replacement of the Newtonian r-modes with hybrids, the structure of the inertial mode spectrum 
in relativistic stars appears to be identical to that in Newtonian stars. That is, there is a one-to-one correspondence 
between the rotationally restored modes of a Newtonian star and those of the corresponding relativistic model. One 
can make a more formal comparison by constructing sequences of relativistic models of increasing compactness, M/R, 
a suitable measure of the importance of relativistic effects. A relativistic star with small M/R agrees well in all of 
its physical properties with the Newtonian model having the same EOS and central density. In each star along the 
relativistic sequence we search for inertial modes using the method described in Appendix We find that for each 
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of the Newtonian modes considered, there exists a family of relativistic modes parametrized by M/R that approaches 
the Newtonian mode as M/R ^ 0. 

To test the accuracy of our numerical code we compared its results with our post-Newtonian solution from Paper II 
in the small M/R regime (see Fig. 0). The solution (II,5.33)-(II,5.40) gives explicitly those axial- led inertial modes of a 
relativistic uniform density star that correspond to the I = m Newtonian r-modes. For these modes the dimensionless 
comoving frequency in the Newtonian star is simply 

(35) 



" {m + iy 

while the post-Newtonian calculation gives Eq. (11,5.33), or, 



(to + 1) 



^ 8(to-1)(2to + 11) fM 



5(2m + l)(2TO-|-5) \R 



(36) 



Fig. shows that the numerically computed eigenvalues agree well with the post-Newtonian solution; the differences 
being of order (A//i?)^ as expected. The numerically computed eigenfunctions also agree well with the post-Newtonian 
eigenf unctions (II,5.34)-(II,5.40). This agreement with our analytic solution gives us confidence that our code is able 
to find the relativistic modes. Thus, we may now explore the fully relativistic regime and, in particular, consider 
modes for which we have not worked out a post-Newtonian solution. Fig. |l| also shows the numerically computed 
eigenvalues for highly relativistic uniform density models. We have used the Newtonian frequency kn to normalize the 
results. This makes it easy to see that all modes clearly approach their Newtonian values as M/R 0. Furthermore, 
the results show that the mode-frequency tends to decrease as the star becomes more compact. It is anticipated 
that general relativity will have this effect The gravitational redshift will tend to decrease the fluid oscillation 
frequencies measured by a distant inertial observer. Also, because these modes are rotationally restored they will 
be affected by the dragging of inertial frames induced by the star's rotation. Specifically, since the Coriolis force is 
determined by the fluid's angular velocity relative to that of the local inertial frame uj{r) = il — uj{r) it decreases, and 
the modes oscillate less rapidly, as the dragging of inertial frames becomes more pronounced. This is, of course, not 
too surprising given ^ 

Next we consider how the relativistic analogues to the Newtonian r-modes are affected by changes in the stiffness 
of the equation of state. We can easily do this by determining the modes for a varying polytropic index, n, at a given 
compactness. Fig. ^ provides such results for M/R — 0.15 which would be a typical value for realistic neutron stars. 
From this data we see that the mode frequencies differ from the Newtonian values by 10 — 15%, depending on the 
stiffness of the equation of state. In order to confirm that these results are typical, we have carried out calculations for 
a much larger set of modes (inluding also polar- led inertial modes). Results for all modes whose Newtonian frequencies 
Kn were listed in Tables 5 and 6 of Paper I, show shifts in the frequency that accord well with the results in Fig. ^. 

We now turn to a discussion of the eigenfunctions. Because the Newtonian studies have found the I = m r-modes 
to be the most unstable to gravitational radiation-reaction, we would like to understand the nature of these modes 
in relativistic stars. As we have already mentioned, these modes become axial-led hybrids in relativistic barotropic 
models. In terms of the spherical harmonic expansion of the fluid displacement, Eq. (|l0|), the statement that the 
Newtonian modes are pure r-modes with I = to means that only the coefhcient Um(r) is nonzero (to lowest order in 
ri). It has the form 



\m+l 

Um{r) = {-) . (37) 



[We have normalized the mode according to Eq. (|3^ ) so that Um{r) = 1 at the surface of the star, r = R] In 
a relativistic barotrope, on the other hand, other coefficients in Eq. (|o|) will be nonzero to lowest order in fi. In 
addition we can examine the nonzero functions hi(r) (= /iq,;) and Hij{r) in Eq. ([T^), the spherical harmonic expansion 
of the perturbed metric. 

Figs. 1^ and ^ show some of these nonvanishing coefficients for the axial-led inertial mode whose Newtonian coun- 
terpart is the I — m — 2 r-mode. The mode is shown for two relativistic models: a uniform density star (n = 
polytrope) and an n = 1 polytrope; both with compactness M / R = 0.15. At this compactness, the mode of the n = 
model has eigenvalue k = 0.5991 and the mode of the n — 1 model has eigenvalue k = 0.5907. [In Newtonian gravity 
the eigenvalue is Kn = 2/(to -I- 1) = 2/3.] Fig. ^ shows the fluid displacement functions Ui{r), Wi{r), and Vi{r) for 
^ < 4 as well as the Newtonian r-mode function U2{r) — (r/R)^ (dashed curve). Meanwhile, Fig. ^ shows the metric 
functions hi{r) and Hi^i{r) for I < 4. The amplitudes of all of these functions are determined by the normalization 
condition, Eq. (Q), and reveal that the relativistic corrections to the Newtonian r-mode eigenfunction are only of the 
order of a few percent even for these highly relativistic models. The fact that h2{r) dominates the perturbed metric 
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reveals directly that this mode couples most strongly to current quadrupole radiation, cf. estimates based on the 
Newtonian mode results |l^, 12 



Given the existence of a large number of inertial modes of a rotating star, it is interesting to ask whether other such 
modes may be driven unstable by gravitational radiation. That this is, indeed, the case was shown in Paper I (see 
also pot). A particular axial-led mode which is interesting since one would a priori expect it to couple strongly to 
current quadrupole radiation is shown in the series of figures |[|^. This mode is the relativistic counterpart of one of 
the Newtonian m = 2 axial-led hybrids studied in Paper I. (In particular, it is the mode whose Newtonian frequency 
is = 0.4669 in the uniform density model and = 0.5173 in the n = 1 polytrope, cf. Paper I.) 

Fig. H shows the Newtonian mode in a uniform density model together with its relativistic counterpart at the nearly 
Newtonian compactness of M/R = 1Q~'^. The Newtonian and weakly relativistic fluid functions Ui{r), Wi{r), and 
Vi{r) for / < 4 arc indistinguishable - the relativistic corrections being of order 10~^, as one would expect. This scale 
for the relativistic corrections is also indicated by the leading metric functions hi{r) and Hi^i{r) and by the fact that 
the functions with I > 4 are smaller than those shown in both panels by a factor of order 10~^ or smaller. 

As discussed in Paper I, it is possible to find exact, analytic expressions for the Newtonian inertial modes of a 
uniform density star. The explicit forms of the fluid functions shown in Fig. ^ are given in Table 3 of Paper I (but 
with a different normalization) . The functions with / > 4 are identically zero for the Newtonian uniform density model, 
so the eigenfunction shown is indeed an exact solution. (For the Newtonian n = \ model the fluid eigenfunctions are 
similar but not identical to those shown and have nonvanishing higher / terms of order 0.5% or smaller.) 

The fact that U2{r) is of order unity relative to the other fluid functions suggests that this mode ought to lead to 
significant current quadrupole radiation. Since it satisfies the CFS instability criterion the mode will be driven unstable 
by gravitational radiation and one might even expect it to make as important a contribution to the gravitational-wave 
driven spin-down of a hot, rapidly rotating neutron star as the much-discussed I ~ m — 2 r-mode. However, a more 
detailed calculation of the growth timescale of the unstable mode does not support these expectations. Instead, one 
finds that the current quadrupole radiation from this mode is negligible compared to higher I multipole moments for 
the Newtonian models. In Paper I, it was shown that the current quadrupole associated with this mode actually 
vanishes identically for the uniform density model and that it is negligibly small for the n = 1 polytrope. The 
relativistic calculation provides us with a way to understand this result since we can now examine the perturbed 
metric in the exterior spacetime for a weakly relativistic model. From the results in Fig. |^ we see directly that the 
metric function h2{r) essentially vanishes outside the star compared to the function h^ir), thus implying that the 
exterior perturbed metric (i.e. the emerging radiation) is dominated by the / — 4 current multipole. 

It is interesting to consider whether this result will still be true for strongly relativistic stars. In Figs. ^ and |^ we 
show the same mode in two stellar models with compactness M/R = 0.15. For this compactness, the mode of the 
n = model (Fig. ^) has eigenvalue n — 0.3879 and the mode of the n = 1 model (Fig. |^) has eigenvalue k ~ 0.4313. 
That is, the frequencies are 17% smaller than the Newtonian ones. For clarity, we display the coefficients of the axial 
and polar-parity terms in the spherical harmonic expansions (u^ and (13) in separate panels. The upper left panel 



of each figure shows the axial-parity fluid functions Ui{r) for Z < 6, while the upper right panel shows the polar parity 
ffiiid fimctions Wi{r), and Vi{r) for I < 5. (Their Newtonian counterparts are also shown for comparison.) The lower 
left and right panels show, respectively, the axial metric functions hi{r) and the polar metric functions Hi^i{r) for 
I <Q. These figures indicate that the relativistic corrections to this mode are only of the order of 1% even for these 
strongly relativistic models. The results also show that the metric function h2{r) is not completely dominated by 
/i4(r) in the exterior spacetime. For both stars, /12 is of order 0.1% at the surface of the star (where, as usual, the 
normalization is fixed relative to U2{R) = 1). This is comparable to /i4 for the n = 1 star and only a factor of 10 
smaller than /14 for the n = star. This is an interesting result since it suggests that the Newtonian prediction that 
the radiation fie ld is dominated by the / — 4 current multipole is not borne out by the fully relativistic calculation. 
Indeed, in Sect. IV C we will see that the I = 2 current multipole dominates over the I = 4 multipole. We will also 
show that the growth timescale of the mode is nevertheless much longer than that of the I = m = 2 r-mode. 

Finally, we present a set of results for the fastest growing unstable polar- led inertial mode in Figs. p[]To|. This mode 
is the relativistic counterpart of one of the Newtonian m — 2 polar-led hybrids studied in Paper I. (In particular, it 
is the mode whose Newtonian frequency is = 1.232 in the uniform density model and = 1.100 in the n ~ 1 
polytrope.) Fig. H shows the Newtonian mode in a uniform density model together with its relativistic counterpart 
at the nearly Newtonian compactness oiM/R— 10"'*. Again, the Newtonian and weakly relativistic fluid functions 
Ui{r), Wi{r), and Vi{r) for I < 3 are indistinguishable — the relativistic corrections and leading metric functions 
hi{r) and Hi^i{r) being of order 10~^. As with the axial-led mode shown in Fig. ^, an exact expression for the mode 
eigenfunction in the Newtonian uniform density star is given in Paper I (Table 4) with the fluid functions vanishing 
identically for / > 3. (Again, for the Newtonian n = 1 model the fluid eigenfunctions are similar but not identical to 
those shown in Fig. ^ and have nonvanishing higher I terms of order 1% or smaller.) 

In Figs. p[]lO| we show the mode in two stellar models with compactness M/R = 0.15. At this compactness, the 
mode of the n — model (Fig. p|) has eigenvalue k = 1.028 and the mode of the n — 1 model (Fig. |l^) has eigenvalue 
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K = 0.9068. This means that the frequencies are again about 17% smaller than the Newtonian results. For clarity, 
we display separately the coefficients of the axial and polar-parity terms in the spherical harmonic expansions ( |lO| ) 
and (^3|). These figures show that the relativistic corrections to this mode are, again, only of the order of 1% even 
for these strongly relativistic models. We can also see that the metric function /i3(r) dominates the perturbed metric 
exterior to the star, which implies that this mode couples most strongly to current octupole radiation. 



IV. RADIATION REACTION TIMESCALES 



Almost all previous estimates of the strength of the r-mode instability are based on Newtonian and post-Newtonian 
calculations. One of the main goals of this work is to determine whether general relativity will have a significant 
effect on the stability of the modes, and, if so, to determine whether it will make them more or less unstable. This 
is obviously a crucial issue. In particular since it is known that relativistic effects change the instability point for 
the f-modes considerably This fact underlines why radiation driven instabilities need to be studied in the full 
framework of general relativity. As far as unstable inertial modes are concerned, there has so far been two studies of 
the associated growth times. Both concern the r-modes of non-barotropic relativistic stars. Yoshida and Futamase 
p9| implemented the near-zone boundary conditions discussed by Lindblom, Mendell and Ipser , while Ruoff and 
Kokkotas |Q solved the complex eigenvalue problem |^ posed by the slow-rotation equations when a condition of 
outgoing radiation at infinity is imposed. These two studies provide useful insights into the growth of the unstable 
r-modes in general relativity, but the problem is still far from well understood. In particular, there are no estimates of 
the radiation reaction timescales for general inertial modes of a fully relativistic star in the current literature. This is 
an unfortunate gap since it seems reasonable to expect that barotropic models (for which all relativistic inertial modes 
are hybrids) will be relevant for rapidly spinning stars. In addition, it was shown in Paper I that, many of the hybrid 
inertial modes are unstable due to the emission of gravitational radiation. Their growth timescales were estimated 
using a post-Newtonian calculation and found to be considerably longer than that associated with the purely axial 
r-mode. It is clearly relevant to ask whether these results remain accurate when relativistic effects are included in the 
analysis. In fact, we have already discussed why this will not be the case for the particular axial- led mode illustrated 
in Figs. 10. 

In this section, we present the first calculation of the gravitational radiation-reaction timescales for general inertial 
modes of fully relativistic barotropic stars. Our method for determining these timescales is significantly different from 
those used in It is close to the post-Newtonian approach in spirit, but based on a fully relativistic analysis. 

When the energy radiated per cycle is small compared to the energy of the mode, E, the imaginary part of the mode 
frequency is accurately approximated by the expression. 

The rate of energy change is determined by a competition between dissipative effects such as viscosity (which tend to 
damp these modes) and gravitational radiation reaction (which drives the unstable modes). To estimate the contri- 
bution to T from gravitational radiation we must calculate both the mode energy and the rate at which gravitational 
waves carry energy away from the star. 



A. Calculation of the mode energy 

We compute the mode energy using the Lagrangian perturbation formalism ^ ^ . The appropriate quantity 
to calculate is the canonical energy, Ec- For canonical displacements Ec agrees with the physical second order change 
in the energy associated with the mode to lowest order in perturbation theory |^] . For a mode with behavior g^(<^t+">-v) 
the canonical energy and angular momentum are related by Ec ~ —{a /m)Jc. Thus, instead of computing the canonical 
energy directly we perform the (much simpler) calculation of Jc using the following expression p6| , p2| , 

Jc^-^^Jric, |(7"^^^i:^e^V^e5 + V'^^X^CphyS - ^e^''^^''ef'P^£^h;s^ph^pj dV (39) 

where * denotes complex conjugation, — e'^V at is the past-directed timelike normal to a t = constant hypersurface, 
dy = is the volume element on that surface and 
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This expression for Jc and the corresponding expression for Ec befow are the canonical angular momentum and 
energy for the physical perturbation, for the real (or imaginary) part of the complex perturbation [ft.Q/3,C"]- Care 
must be taken in comparing their values to those of other papers ^ that define Jc (or Ec) for a complex 

perturbation to be the sum of its values for the real and imaginary parts. 

For our slowly rotating star, we need to compute the mode energy only to leading order in f2. Taking into account 
the ordering (|l7|), a laborious but straightforward calculation gives. 



Ec 



m 



1 

r 



dV. 



(40) 



[Note that it is necessary to use the fact that ^" is a canonical displacement to reduce expression ( p}S| ) to this form.] 
We now substitute for ^" and ha(3 their spherical harmonic expansions (^) and ( |l3|) and perform the angular 
integration. This leaves us with an expression for the canonical energy involving the variables Ui{r), Vi{r), Wi{r) and 
hi{r) associated with the hybrid mode eigenfunctions. 



with e the average energy density and 



E 

l—m 



Ec 



_ ^r{y' + A') 



am 



Wf + l{l + l)Vf + /(/ + l)Ui{Ui +hi)}d 



R 



(41) 



(42) 



Given one of our numerical solutions for a particular eigenmode it is straightforward to perform this radial integral 
numerically and compute the mode energy. It is also straightforward to show that (^) reduces to Eq. (64) of Paper I 
in the Newtonian limit, once we account for the difference between the inertial and the rotating frames and for the 
fact that (1,64) involves the sum of the real and imaginary parts of a complex perturbation, giving it an extra factor 
of two relative to (^T]). 

B. Calculation of the radiated power 



We now compute the rate at which energy is emitted in gravitational waves, simplifying the analysis by using a 
gauge- invariant expression (derived in Appendix ^) for the asymptotic power radiated by an axial mode. We otherwise 
closely follow a calculation of Ipser ||50|| (our bibliography notes two minor corrections to his equations.) 

Recall that the perturbation equations we derived in Sect. II B are relevant only in the near zone. Because the 



mode frequency is proportional to the star's angular velocity, however, the slow rotation approximation allows one to 
extend the near zone arbitrarily far beyond the star by making the angular velocity sufficiently small (see Table |). 
Hence, to find the more general equations valid in the wave zone as well as in the near zone, we need only concern 
ourselves with metric perturbations of the exterior vacuum spacetime. Furthermore, the equations relevant here are 
those governing the perturbations of a spherical background. The leading rotational terms are of order wr, and outside 
the star the metric function uj{r) has the form, 



uj[r) 



2J 



(43) 



with J the angular momentum of the star p4| . Such terms can be neglected in the slow rotation approximation 
because 



r I \ r 



(44) 



Hence, the relevant components of the perturbed Einstein equation [SG^f = and 6Gg — 0), are 



K,i ^ 



2M 



hi,. 



, 2A/\ , , / 2A/\ 2M , 
= 1 /ii i + 1 hi^i - iaho,i, 



(45) 
(46) 
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where we have restored the "0" subscript on the metric function /io,i(= ^i) to distinguish it from hij. 

These equations may be combined to give wave equations for either Hq i or hi i. It is not difficult to show that the 
resulting equation for ho^i reduces in the near zone to Eq. (p^, which we used to impose boundary conditions on our 
interior solution. The wave equation for hi^i is, 



r \ i~ J \ r J ^'^ 



(P + l- 2) 



2M\ AM / 5M 



r ) 



(47) 



Following Ipser |5^, we make use of the fact that the near zone extends well into the nonrelativistic region (see 
Table |): Eq. ( ^7| ) is easily solved in the nonrelativistic limit (M/r — > 0) and the solution will be valid both in the 
wave zone and in the outer part of the near zone. Thus, we will be able to impose outgoing-wave boundary conditions 
at r — > oo, and still match to our near-zone solution, Eq. (p9[). In the nonrelativistic zone, Eq. ( ^ ) is simply 



[1^ + 1-2) 



with solution 



hi,i = 



(/-1)(2;-1)!! 



hi I = 



[Ani{ar) + Bji{ar)] 



(48) 



(49) 



where A and B are constants, n/ and ji are spherical Bessel functions and the overall normalization has been chosen 
for later convenience. This solution for hi^i gives rise to a solution for /iq,; using Eq. (^6|), which, in the nonrelativistic 
zone becomes simply, 



icr/io,i = h'li. 

In the outer part of the near zone (or << 1) we have 

B (/ + 2)((Tr)2'+i 



hoA ~ A ( - I <! 1 



A(/-1)(2;h-1)[(2Z-1)!!] 



[1 + 0(^^2)] 



In the wave zone (or >> I) we have 

(ai?)'+i 



(/- 1)(2; - 1)!! \R 



( — ) [Asin((Tr - Itt/2) + B cos{ar - lTr/2)] . 



(50) 



(51) 



(52) 



On this wave zone solution we now impose the boundary condition that the radiation be purely outgoing. With 
time dependence e"^* this requires that B = iA. Our solution then becomes. 



wave zone: 



near zone: 



iA{aR) 



i+i 



{i-i){2i-iy.\ \RJ 



41 e-t" 



-lir/2) 



(53) 



A 



R 



z(? + 2)(f7r)2'+i 



^|[l + 0(aV)]. (54) 



(/-1)(2Z + 1)[(2;- 

The normalization constant A is fixed by matching Eq. (p4) in the outer (nonrelativistic) part of the near zone to the 



solution (|29|), which is valid throughout the near zone. This gives A = hi^o and justifies our neglect in Sect. II B of the 
singular solution to Eq. (|2^ ) by the fact that the singular solution is indeed of order (crr)^'+^ relative to the regular 
solution. The constant hi ^ is, in turn, fixed by the matching to the interior solution, Eq. (^2|), whose normalization 
ultimately is set by Eq. ([3^). 

We now have a solution for our axial-parity metric perturbation valid in the entire domain r e [0, oo). To compute 
the rate of energy radiation we use the gauge-invariant expression 



^) = -7^ii-m + m + ^) lini 

at / SZTT r^co 



(55) 
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where 



fco ^0 - -j^dth^. (56) 



In Regge- Wheeler gauge we have /12 — 0, ko — ho, and hence we obtain 

dE ^ l{l + l){l + 2){aRf+^\hi^o\^ 
dt ~ 327r(/- 1)[(2Z- l)!!]2i?2 ' 

C. The gravitational radiation reaction timescale 



(57) 



We now combine our expressions ( ^l|) for the mode energy and (|57| ) for the power radiated to find the gravitational 
radiation reaction timescale (|38|) . Restoring factors of G and c, we write the timescale as 



^'^'^ l>2 



Ton f::iTi \7rGeJ 

c /3G 

with the integral / defined by Eq. ([42[). 

Having obtained this expression, we are in a position where we can evaluate the growth timescales associated with 
the unstable inertial modes of relativistic stars. As a quick check on our calculation, we plug into Eqs. ( ^2| ) and 
( p9| ) our post-Newtonian solution (II,5.33)-(II,5.40) for the corrections to the / = m Newtonian r-modes of a uniform 
density star. The resulting timescale associated with the I — ni — 2 multipole agrees with previously published results 

m 

1.4Mf^W R V fnGe 



where 



1 ^ c /3GMY ?(/ + !)(/ + 2)^(^-m)2'+i|fe;^o|2 
R\4cmj 16(/- 1)[(2;- 1)!!]2 / 



n 



ro. = 1.56s . (60) 

(See, for example. Table 8 of Paper I.) 

As a further check on the calculation we verify that, for weakly relativistic models, the timescale (^8|) exhibits the 
same scaling with the stellar parameters as the r-mode growth timescale estimated from the Newtonian calculations. 
The Newtonian timescales for the I = ni r-modes exhibit the following scaling with the neutron star mass and radius 
011: 

GM fGM\-^'+'^ fnGeY^''> 
= j j (61) 

where T; is a constant that depends only on I and the star's EOS. with To test our formula for the growth timescale 
we keep the baryon mass, Mb, fixed at 1.4Mo and set fl^ — nGe. One would then expect a log- log plot to clearly 
reveal that the timescale depends on the star's compactness as (Af/i?)~'^'+^^ for low M/R. That this is, indeed, the 
case can be seen from the data in Fig. Ill], which compares the Newtonian and relativistic growth timescales of the 
modes whose Newtonian analogues are the first five I — m r-modes. Fig. |l^ illustrates our results in a different way 
by showing the dependence on the polytropic index, n, of the timescale for the I = m = 2 r-mode and its relativistic 
counterpart. 

These two figures suggest that for highly relativistic stars, the relativistic calculation tends to give a slightly longer 
growth timescale than that of a Newtonian star with the same EOS, baryon mass and compactness. This suggests that 
general relativity tends to stabilize the modes making the r-mode instability slightly weaker than previously expected 
(in accordance with the results for nonbarotropic stars |3^). That this should be the case is natural: All inertial 
modes have relatively low frequencies since a fl. This means that the associated gravitational waves will suffer 
significant backscattering by the spacetime curvature as they escape to infinity. As the star becomes increasingly 
compact more of the curvature potential, which can be approximated by 

V.fl-'-M.) (^-^)(^ + ^) (62) 
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Model 


M/R 


M [Mq] 


R [km] 


Pc [10^5 g/cm=*] 


ec [ICP^ erg/cm^] 


Pc [lO^"* dyne/cm^] 


same 


NO 


0.1650 


1.400 


12.53 


0.3380 


3.038 


2.506 


- 


1 


0.1440 


1.269 


13.02 


0.2734 


2.457 


2.506 


Pc 


2 


0.1490 


1.264 


12.53 


0.3053 


2.744 


2.940 


R 


3 


0.1538 


1.260 


12.10 


0.3380 


3.038 


3.408 


pc-, Ec 


4 


0.1650 


1.248 


11.17 


0.4251 


3.820 


4.763 


M/R 


Nl 


0.1650 


1.400 


12.53 


1.1120 


9.994 


8.245 


- 


5 


0.1390 


1.306 


13.87 


0.7819 


7.852 


8.245 


Pc 


6 


0.1493 


1.300 


12.85 


0.9813 


9.994 


11.746 




7 


0.1529 


1.297 


12.53 


1.0601 


10.855 


13.272 


R 


8 


0.1552 


1.296 


12.33 


1.1120 


11.426 


14.319 


Pc 


9 


0.1650 


1.291 


11.55 


1.3552 


14.150 


19.702 


M/R 



TABLE II: Comparison between i) a fiducial Newtonian uniform density star (labelled NO) and four relativistic uniform density 
models (labelled 1 to 4), and ii) a fiducial Newtonian n = 1 model (Nl) with five relativistic models with the same poytropic 
index (5-9). All of the relativistic models are constructed so as to have a baryon mass of 1.4M0 and also to agree with the 
fiducial Newtonian model on one of their other physical characteristics (indicated in the final column). Gravitational radiation 
reaction timescales for unstable modes of these stellar models are presented in Tables Fl| and IV 



cf. Eq. (^3)> is unveiled. Hence, one would expect low- frequency modes of oscillation to radiate less efficiently as the 
star becomes more compact. 

We should point out that it is somewhat misleading to do a direct comparison of Newtonian and relativistic models, 
because there is no one-to-one correspondence between the two. Although a weakly relativistic polytrope agrees well 
in all of its physical characteristics (mass, radius, compactness etc.) with the Newtonian model having the same 
polytropic index and central density, this is not true for strongly relativistic models. Newtonian and relativistic stars 
with the same polytropic index can be constructed so as to agree with respect to two of their physical characteristics, 
but in general they will differ on the others and this will affect the growth timescales of unstable modes. 

We would nevertheless like to quantify how our fully relativistic radiation timescales differ from the Newtonian ones. 
To account for the lack of one-to-one correspondence between Newtonian and relativistic models, we have constructed 
a number of different relativistic stars that agree in two of their physical characteristics with a given Newtonian 
polytrope of the same index, n. The fiducial Newtonian models are chosen in such a way that a 1.4Mq star has a 
radius of 12.53km. This facilitates comparison with results from the literature, which made use of stars with these 
parameters ||ll|, |l2|, ^ It is natural to choose one of the physical characteristics on which the relativistic 

and Newtonian stars agree to be the baryon mass, Mb, (i.e. the rest mass) of the star. (For relativistic stars the 
baryon mass is slightly higher than the total gravitational mass, whereas these quantities are the same for Newtonian 
stars.) For the other physical characteristic on which the stars agree we choose such quantities as the radius, central 
pressure, central energy density and so on. We list the characteristics of these various models in Table ^ In the 
table we compare a number of relativistic uniform density (n = 0) stars with our fiducial Newtonian uniform density 
star, labelled NO. Similiarly, we compare a set of relativistic n = 1 polytropes with our fiducial Newtonian n = 1 
polytrope, labelled Nl. All of these stars are constructed so as to have a baryon mass of 1.4Mq and also to agree 
with the fiducial Newtonian model on one of their other physical characteristics. The models are listed in order of 
increasing compactness and the final column of the table indicates the second physical characteristic on which the 
stars agree. 

The effect of this ambiguity in comparing the relativistic and Newtonian models is indicated in Table [II, which 
presents the gravitational radiation reaction timescales for the fastest growing I = m Newtonian r-modes and their 
relativistic hybrid counterparts. (Some of these Newtonian timescales have been computed in previous work 
[l^ , p8| p9| , p4|.) We see that, depending on which physical characteristics one chooses to equate, the timescales for 
the analogous mode of the Newtonian and relativistic models can differ by as much as an order of magnitude, with 
the relativistic mode generally having the longer growth time. The fastest relativistic growth times are obtained by 
equating the compactness, M/R, of the relativistic and Newtonian models. In this case, the relativistic growth times 
are typically weaker than the corresponding Newtonian growth times by only a factor of a few. 
In connection with the results listed in Table 



III 



it is relevant to make two observations. First of all, the tabulated 
data suggest that for each / there is a variation of about a factor of two in f; between the various relativistic models (for 
the same equation of state). It is relevant to point out that had we instead tabulated the combination f;A/'+'^/i?'+'^ 
(i.e. accounted for the expected scaling with mass and radius) then the variation between the models would have 
been much smaller. For I = 2 we would have found a variation of about 6% between the uniform density models. 
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Model M/R 


I = m = 2 I — m — 3 I — m = 4 I = m = 5 I — m — 6 


NO 0.1650 


-1.56 X 10" -1.17 X 10^ -8.79 x 10^ -6.19 x 10^ -4.11 x 10=^ 


1 0.1440 

2 0.1490 

3 0.1538 

4 0.1650 


-3.84 X 10" -5.60 X 10^ -7.43 x 10^ -8.87 x 10^ -9.76 x lO" 
-3.27 X 10" -4.70 X 10^ -6.14 x 10^ -7.19 x 10^ -7.76 x W 
-2.82 X 10" -4.01 x 10^ -5.16 x 10^ -5.94 x 10^ -6.31 x W 
-2.03 X 10" -2.82 X 10^ -3.52 x 10^ -3.92 x 10^ -4.01 x W 


Nl 0.1650 


-3.26 X lO" -3.11 X 10^ -2.84 x 10^ -2.37 x 10^ -1.81 x lO" 


5 0.1390 

6 0.1493 

7 0.1529 

8 0.1552 

9 0.1650 


-1.09 X 10^ -2.23 X 10^ -3.80 x 10^ -5.57 x 10'' -7.32 x lO'^ 
-7.95 X 10" -1.59 X 10^ -2.63 x 10^ -3.72 x 10'' -4.70 x 10^ 
-7.15 X 10" -1.42 X 10^ -2.33 x 10^ -3.25 x lO" -4.06 x lO'^ 
-6.70 X lO" -1.33 X 10^ -2.17 x 10^ -3.01 x lO" -3.71 x 10^ 
-5.15 X 10" -1.01 X 10^ -1.61 x 10^ -2.18 x lO" -2.61 x 10^ 



TABLE III: Gravitational radiation reaction timescales, fj, in seconds for unstable (f; < 0) modes of the stellar models listed 
in Table The growth timescales listed here are those of the fastest growing I — m Newtonian r-modes and their relativistic 
counterparts. In general, the relativistic models (1-9) produce longer timescales than those of their Newtonian analogues (NO 
and Nl), suggesting that general relativity tends to stabilize the modes slightly. 



while the result for the polytropes vary by about 11%. This indicates that the variation in the growth timescale 
between models 1-4 and models 5-9 is mainly due to the differences in mass and radius. The second feature worth 
noticing from the data in Tabic is that, while the relativistic timescales differ from the Newtonian ones by only a 
factor of 2-3 for the quadrupole mode, the difference increases with I. For example, for / = 6 the difference is at least 
an order of magnitude. This result can likely be explained in terms of backscattering from the curvature potential 
in the exterior spacetime. From ( |6^ we see that the "height" of the potential increases as {I — 1)(^ + 2). Thus we 
would expect the difference between our fully relativistic results and the Newtonian ones to increase with / roughly 
as {I — 1)(^ -I- 2)/4 (after rescaling with the result for the quadrupole mode). Our numerical results are in reasonable 
agreement with this expectation. 

Finally, we want to confirm the expectation that the analogue of the I = m = 2 Newtonian r-mode remains the 
fastest growing unstable mode also when the growth times are estimated in full general relativity. That this is the case 
can be seen from Table IV where we list the growth timescales for a number of unstable inertial modes. We compare 
the timescales of modes of the Newtonian models (NO and Nl) with those of the corresponding modes of relativistic 
models 4 and 9. (These are the models with the same compactness as their Newtonian counterparts, and which lead 
to the fastest growth times, cf. Table HI). For each mode considered, we list (in four consecutive rows) the data 
for the four different stellar models — all of which have baryon mass Mb = 1.4M0 and compactness M/R = 0.165. 
We list the frequency of the mode in each star as well as the growth timescales associated with the various current 
multipole moments of the mode. We also list (enclosed in parentheses) the timescales associated with some of the 
inertial mode mass multipoles of the Newtonian n = 1 polytrope. These mass multipoles are higher order in than 
can be computed within our slow rotation formalism. However, they have been calculated by Yoshida and Lee using 
a self-consistent third order 

Newtonian formalism (see Table 4 of Ref . ) and we list them here simply for ease of comparison with our new 
results. The point is to compare the new relativistic timescales with the previously published Newtonian timescales 
(see, in particular. Tables 7-9 of Paper I, Table 4 of |||| and Table 1 of @). 

As with the timescales listed in Table HI , most of the relativistic growth times presented in Table |^ are basically 
unchanged compared with their Newtonian analogs. There are, however, cases where there is a dramatic difference 
between the Newtonian and the relativistic results. The best example of this is provided by the fourth mode listed 
in Table This is the axial-led inertial mode presented in Figs. in Sect. Ill, whose current quadrupole moment 
vanishes (or nearly vanishes) in the Newtonian models. In Paper I, it was argued based on the Newtonian slow-rotation 
calculation that the growth of this mode is dominated by its Z = 4 current multipole. However, by including rotational 
corrections to higher order in O, Yoshida and Lee | p9[ were able to compute the I = 3 mass multipole and found that 
it drives the mode on an even shorter timescale. Now, with the inclusion of relativistic corrections to the mode, we 
see that the growth of the mode is, in fact, dominated by the I = 2 current multipole — as one would have expected 
in the first place. This mode is significantly more unstable in general relativity than in the Newtonian calculations. 
However, it's growth time is nevertheless still much longer than that of the mode whose Newtonian analogue is the 
I = m = 2 r-mode. 



To conclude: The data presented in Tables III and IV suggest that the relativistic corrections to the timescales 



are not large enough to alter the standard picture of the gravitational-wave driven instability of sufficiently hot and 
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m 


Parity Model 


K 


T2 






T4 


^5 


-1 
i 


a 


NO 


n 1 on 
U.Di/U 






n Id 1 a6 

— y. (a X iu 




— OO 






4 


0.5008 






-5.21 X 10*^ 




-4.19 X 10^" 






Nl 


0.6906 ( 


-2.46 X 10^) 




-1.25 X 10** 




-1.22 X 10^" 






9 


0.5630 






-4.65 X 10^ 




-2.06 X 10^^ 


o 
I 


a 


NO 


0.666/ 


— i.OD X iO 






— OO 








4 


0.5903 


-2.03 X lO" 






-2.09 X 10® 








Nl 


0.6667 


-3.26 X 10" 


( 


-3.49 X 10^) 


— oo 








9 


0.5796 


-5.15 X 10" 






-8.51 X 10* 






p 


NO 


1.2oi9 






— 4. / / X iO 




— OO 






4 


1.0039 






-2.41 X 10* 




-6.07 X 10^^ 






Nl 


1.1000 ( 


-1.71 X 10^) 




-3.37 X 10* 




-3.13 X 10" 






9 


0.8780 






-3.41 X 10* 




-1.56 X 10^2 




a 


NO 


0.4069 


— oo 






Q OO 1 A^ 

— 3.00 X iO 








4 


0.3786 


-1.33 X 10* 






-1.16 X 10" 








Nl 


0.5173 


< -10i« 


( 


-8.39 X 10*) 


-1.85 X lO'' 








9 


0.4206 


-8.29 X 10^ 






-7.47 X 10" 




6 


a 


NO 


A C AAA 

O.oOOO 






1 1 ^7 1 a1 

— i.i r X iO 




— oo 






4 


0.4278 






-2.82 X 10^ 




-1.72 X 10® 






Nl 


0.5000 






-3.11 X 10^ 


(-1.88 X 10^) 


— oo 






9 


0.4259 






-1.01 X 10^ 




-1.04 X 10® 




p 


NO 


1 ACOO 








O AA 1 a4 

— z.OO X iO 








4 


0.8438 








-3.93 X 10* 








Nl 


0.9049 




( 


-8.62 X 10^) 


-2.71 X 10* 








9 


0.7213 








-9.69 X 10* 






a 


NO 


0.3779 






— oo 




-7.67 X 10^5 






4 


0.3057 






-7.20 X 10* 




-4.31 X 10" 






Nl 


0.4126 






< -10^° 


(-5.30 X 10^) 


-3.97 X 10" 






9 


0.3369 






-1.03 X 10* 




-3.36 X 10^ 



TABLE IV: Gravitational radiation reaction timescales in seconds for unstable (ji < 0) rotational modes of Newtonian and 
relativistic stellar models. For each mode, we compare data from four different stellar models: a Newtonian uniform density 
star (model NO), a relativistic uniform density star (model 4) a Newtonian n = 1 polytrope (model Nl) and a relativistic 
n = 1 polytrope (model 9). All of these models have baryon mass Mb ~ 1.4M0 and compactness M/R = 0.165. (Hence the 
Newtonian models agree with the canonical model typically used in the literature with mass 1.4M0 and radius 12.53 km.) We 
list the azimuthal index, m, the dimensionless comoving frequency, k, and the parity of the mode (i.e., whether it is an axial-led 
or a polar-led hybrid) as well as the current multipole radiation timescales computed to lowest order in our slow-rotation 
formalism. For convenience, we also show (in parentheses) the mass multipole radiation timescales computed by Yoshida and 
Lee for modes of the Newtonian n — 1 polytrope. (These mass multipoles, which we are unable to compute within our 
slow- rotation formalism, were computed using a self-consistent third order calculation.) 

rapidly rotating neutron stars. The fastest growing mode (by at least an order of magnitude) is the axial-led inertial 
mode corresponding to the I = m ~ 2 Newtonian r-mode. And although the actual growth timescale of this mode 
is uncertain due to the uncertainty in the neutron star EOS, it is unlikely to be significantly shorter than has been 
estimated here. 
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V. CONCLUDING REMARKS 



In this paper we have studied the inertial modes of rotating relativistic stars. Numerical results were presented for 
the mode-eigenfrequencies and the associated eigenfunctions of barotropic models. These results were shown to be 
in good agreement with results in the literature (in particular in the post-Newtonian limit) , and provide a significant 
improvement on previous studies as far as the strongly relativistic regime is concerned. 

We also analyzed the rate at which these modes radiate gravitationally. In particular, we studied the growth 
timescale of various modes that are unstable due to the emission of gravitational waves |^ . Our calculation was 
based on two ingredients: The energy associated with the mode oscillation was determined as the canonical energy 
defined by Friedman |^ , while the gravitational- wave luminosity followed from an analysis parallel to that of Ipser 
|^o| . By combining these two quantities we arrived at the required damping/growth timescale due to gravitational- 
wave emission. Our approach to the problem is novel, and differs from the methods previously used to estimate the 
corresponding timescales for the modes of non-barotropic stars [^8[ ^ . Nevertheless, the final results show the same 
qualitative behaviour. Most notably, the post-Newtonian estimates for the growth rate of an unstable inertial mode 
are found to be surprisingly accurate even for strongly elativistic models. Still, the results show a deviation from the 
post-Newtonian estimates as the star reaches compactness similar to that expected of a neutron star M/R ~ 0.1. 
Then the efficiency of radiation reaction tends to decrease. This result is likely due to the fact that the low-frequency 
waves from an inertial mode experience enhanced backscattering by the spacetime curvature as the star becomes 
increasingly compact (recall that the spacetime of a spherical star has a curvature potential barrier with a peak in the 
region R/M ^ 3). This means that general relativistic effects tend to stabilize the inertial modes. This is in contrast 
to the results for the instability associated with the acoustic f-modes. As was shown by Stergioulas and Friedman [Q, 
the f-modes are significantly destabilized by relativistic effects. 

Despite significant progress in the last few years, it is still not clear to what extent the gravitational-wave driven 
instabilities in rotating compact stars are of astrophysical relevance. In particular, we do not yet have a clear answer 
to the question of whether the unstable modes may lead to detectable graviational waves or whether they limit the 
spin of nascent neutron stars or of old neutron stars spun up by accretion. However, it is important to realize that 
serious theoretical challenges need to be overcome if we want to make further progress in this area of research. For the 
unstable r-modes, key questions include the role of complex (not well understood) interior physics, e.g. the strength 
of hyperon bulk viscosity p^, |23 , and the saturation amplitude set by nonlinear coupling ||2l| . 



APPENDIX A: THE SPECTRAL METHOD USED TO SOLVE THE EIGENVALUE PROBLEM 

To solve the set of equations (|l^)-(|2^) for the inertial modes of a barotropic relativistic star, we use a variant of 
the method developed in Paper I for the analogous Newtonian hybrid/inertial mode problem p^ . We express our 
equilibrium and perturbation variables as a sum over a set of basis functions and substitute these series into our system 
of differential equations. This results in a system of algebraic equations for the expansion coefficients, which may 
then be solved using standard linear algebra techniques. In Paper I, we made use of power series expansions and were 
able to accurately compute the Newtonian hybrid mode eigenvalues and eigenfunctions. However, for the relativistic 
problem we have found it necessary to use a basis of orthogonal functions and we work instead with Chebychev 
polynomials (or, more specifically, "Type I" Chebychev polynomials ]55[| ). In other words, we use a spectral method 
to approach the problem. 

Since the Chebychev polynomials are naturally defined on the domain [— 1 . 1] , we define a new coordinate 

y = 2Q-l (Al) 

which maps the interior of the star, < r < i?, to the required domain. The Chebychev polynomial of degree i is 
then given by 

Ti(j/) = cos(i arccosy). (A2) 

We begin by expressing our equilibrium variables in terms of this Chebychev basis. We construct an equilibrium 
star using standard integration recipes and then use Chebychev approximation |Q to find a Chebychev series that 
accurately fits each of our equilibrium variables. In other words, we take all of the background variables appearing in 
Eqs. (|l9|)-(|2^), such as e(r), p{r), uj{r) etc., and represent each of them by a Chebychev series of the form. 



i=0 



(A3) 
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where the coefficents, 6^, are determined from the Chebychev approximation algorithm (see Numerical Recipes p6[). 

Having constructed Chebychev series for our (known) background variables, we write each of our (unknown) per- 
turbation variables in terms of an expansion in Chebychev polynomials of the form, 



Fiir) 



r \ '+9 

'r. 



^fi,i Ti{y) - ^fi^o 

.1=0 



(A4) 



The factor {r/Kf^'^ provides for the condition of regularity at the origin, which requires the perturbation variables 
to vanish as an appropriate power of r as r — > 0. The axial parity variables Ui{r) and hi{r) have q = \, while the 
polar parity variables Wi{r) and Vi{r) have q — Q. 



If we substitute the various Chebychev series represented by (A_3) and (A4) into our perturbation equations, each 
term in these equations will take the form of a product of two Chebychev series. That is, the generic term in our 
perturbation equations will have the form B{r)Fi{r) with B{r) a known background function and Fi{r) an unknown 
perturbation variable. We would like to be able to write such a product as a new expansion in Chebychev polynomials. 



This can be accomplished using the identity 2TiTj — Ti^j + T\i-j\, which follows from Eq. (A2) and standard cosine 



identities. After some careful rearrangement of terms we find for the product of the series (^ and (0) 



B{r)Fi{r) 



r \ '+9 

'r. 



.i=0 



where 



with 



e(fc) 



bi+j + B(j - l)fo|j-j| 



for fc < 
for fc > 



(A5) 



(A6) 



(A7) 



We also need an expression for the derivatives of our perturbation variables in terms of the Chebychev expansions 



(A4). If we define the Chebychev expansion for the derivative of F;(r) as follows: 



R 



dr 



\R 



-(i+q) 



Fi 



i=0 



(A8) 



and then make use of standard identities involving Chebychev polynomials 
coefficients fi,i of this series are related to the coefficients fi,i of (A4) by 



it is not too difficult to show that the 



fi,i - fi,i+2 = 4(i + l)/;,i+i. 
Our method of solution is as follows: We expand all of the quantities appearing in Eqs. 



(A9) 



J|)-(|2^) in Chebychev 

series, and substitute these exp ansions into the equations and into the boundary and matching conditions. We then 
use the formulas (A5) and (Ra) to express the resulting equations as a linear algebraic system of the form 



Ax = 



(AlO) 



where A is a known matrix that depends nonlinearly on the parameter k and x is a vector whose components are the 
unknown coefficien ts in the Chebychev series for the variables hi, Ui, VJ, Wi and their derivatives. 

To satisfy Eq. ( AlO ) we must search for those values of k for which the matrix A is singular; that is, we must 
find the zeroes of the determinant of A(k). Since A is infinite dimensional, we must truncate our spherical harmonic 
expansions ( [To| ) and ( p^ ) at some maximum index l^^^ and also truncate our Chebychev expansions (A3), (A4) and 
(A8) at some maximum index i,-^^^. The resulting finite matrix is band diagonal. To find its zeroes we use standard 
linear algebra and root finding routines. We then check for convergence of these eigenvalues as we increase l^^^ and 

«max- 

The eigenfunctions associated with these eigenvalues are determined by the perturbation equations only up to 
normalizati on. G iven a particular eigenvalue, we find its associated eigenfunction by replacing one of the equations in 
the system (AlO) with the normalization condition (|3^). Since this eliminates one of the rows of the singular matrix 
A in favor of the normalization equation, the result is an algebraic system of the form 



Ax = 6, 



(All) 



where A is now a non-singular matrix and 6 is a known column vector. We solve this system for the vector x and 
reconstruct the various Chebychev expansions from this solution vector of coefficients. 
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APPENDIX B: GAUGE-INVARIANT EXPRESSION FOR THE LUMINOSITY 



We present here a brief derivation of the gauge-invariant expression (55) for the rate at which energy is radiated to 
future null infinity by an axial mode. One first writes the energy radiated in an asymptotically regular gauge 
and then observes that the gauge-dependent quantity in the expression can be replaced by a quantity that is gauge 
invariant. 

We begin with standard expressions for dE/dt. Written in terms of the leading term ct° in the asymptotic shear 
of outgoing null geodesies, in outgoing null coordinates {u, r, 9, ip) of a flat asymptotic metric f7a/3, the instantaneous 
power radiated is 

^ ^ \dta\u,e,^f dn. 



dt 



An 



(Bl) 
1, the 



where / :— lim / , with Sr a sphere of constant r and u. In an asymptotically regular gauge p 

/ r — >oo I c — ' ' — 

components of the metric perturbation in an orthonormal basis {e^^)} fall off like l/r or faster in an outgoing 

null direction, and (for a real perturbation), Eq. (Bl) takes the form. 



dE 
'dt 



1 

16^ 



{dth 



1 



{dth 



dth 



dn. 



(B2) 



Because this expression is identical to the flat-space expression in terms of the deviation h^f} of the metric from 
its Minkowski value, the radiation field can be treated as a perturbation of flat space. In particular, the parts 
of ha[} belonging to different representations of the rotation group decouple, and we can restrict consideration to 
perturbations hafs belonging to an {l,m) representation with axial symmetry. The gauge invariance of the vector 
in (55) arises from this decoupling. 

In Regge- Wheeler notation Goj, a (complex) axial {l,m) perturbation has the form 



htA = -/lo* 



hrA = — 



hAB = h2 XlAB 



(B3) 



with A, B indices on the sphere. Here, as in Eq. (p^, the axial vector harmonics have the form ^Ya — ^a^bYi"^, with 
^ji-^, = — sin 9, while the axial tensor harmonics have components 



-1 



1 



Xe9 



Xe^ = -isin0(a2 - coted^ - -^X^Yi" 



smy- 2"--'"'"" sin' 9"'^''' ' ^^^^ 

A gauge transformation associated with a vector field changes hap by \7aCl3 + V/jCa. The transformation 
preserves the axial {l,m) representation to which the physical perturbation belongs if and only if (a is itself an axial 
vector on the sphere: (t — Cr — 0, 



Ca = c*: 



Then 
and 



htA htA + dtC 



lA^ 



T-rA 



I A 



hrA+r'dr{^0 '^Ta, 



(B5) 



hAi 



ho^ ho- dtC, 



h,^h,-r'dr{\0, 



hAB + CxTab, 
h2 + 2C, 



(B6) 



implying |57 



that the vector /c^, whose nonzero components are 



fco = h.Q + i(9fft.2, 



ki=h^ + lr'dr{r-'h2), 



(B7) 



is gauge invariant. 

When an axial perturbation is written in an asymptotically regular gauge, expression (B2) is valid. For a complex 
perturbation, 



dE , 



(h) 



dE , 
'di 
1 

16^ 
1 



dE 
'dt 



\dth2 



1 I |2 

T \X(e){e) - X{v){v)\ 



dn 



- — {1-1)1(1 + + 2) lim 

lOTT r^oo 



dth2 



2r 



(B8) 



21 



Finally, noting that, in this gauge, /12 = 0{r), Hq — 0(1), we have 



2r r 



^ih)^^-^{i~mi + m + 2) lim 

at iDTT r^oo 



(B9) 



for a complex perturbation. A real axial mode then radiates average power 



^m) = (^{^h)\ mi + m + 2) ^ijm 



dt 



dt 



(BIO) 



with the limit, as above, taken along a radially outgoing null direction. Although, for an unstable mode, the perturbed 
metric and its gauge invariant quantities blow up exponentially at spatial infinity^ the power radiated to null infinity 
is well-defined and finite. 
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Correcting an earlier calculation by Andersson |H] that was inconsistent in that it included some terms of order Q,^, while 
neglecting others. 
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FIG. 1: Left panel: Shift in the normalized frequencies, k/kn, of the axial-led hybrid modes whose Newtonian counterparts 
are the / = m = 2, 3, 4 and 5 r-modes. These frequencies have been computed for a uniform density model (n = polytrope) 
in the small M/R regime. The frequencies computed using our fully relativistic code (solid curves) agree well with the post- 
Newtonian solution (dashed lines) and deviate, as expected, by a correction of order (M/R)^. Right panel: The corresponding 
mode frequencies in the strongly relativistic (large M/R) regime. 
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FIG. 2: Shift in the frequencies of the modes shown in Fig. ^ for a sequence of stars of varying polytropic index, n. Each model 
in the sequence is chosen to have the same compactness, MjR — 0.15. 




FIG. 3: The relativistic axial-led inertial mode (solid curves) whose Newtonian counterpart is the I = m = 2 r-mode (dashed 
curve). We show the fluid functions Uiir), Vi{r) and Wiir) for / < 4 in a uniform density model (left panel) and an n = 1 
polytrope (right panel); each with compactness M/R = 0.15. The vertical scale is set by the normalization of U2{R) = 1 [cf. 
Eq. (64)]; however, the other functions have been scaled up by a factor of 10 to make them visible in the figure. The fluid 
functions with I > 4 are of order 0.05% or smaller and are not shown. 
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FIG. 4: The metric functions hi{r) and Hi^i{r) for I < 4 for the same mode as in Fig |3[ The vertical scale is again set by the 
overall normalization of U2{R) = 1, as in Fig ^. The metric functions with I — 3 and 4 have been scaled up to make them 
visible in the figure, while those with I > 4 are of order 0.005% or smaller and are not shown. 




FIG. 5: An axial-led hybrid mode of a relativistic uniform density barotrope (solid curves). Also shown (dashed curves) is the 
Newtonian counterpart to this mode: the m — 2 axial- led hybrid with frequency kn ~ 0.4669 (see Paper I). The Newtonian 
and relativistic fluid functions are indistinguishable because the mode is shown in the weakly relativistic regime, for a star with 
compactness M/R = 10"*. The left panel shows the fluid functions Ui{r), Vi(r) and Wi{r) while the right panel shows the 
metric functions hi{r) and Hi^i{r), all for I < 4. All of the functions are shown to scale; thus the scale of the right panel reveals 
the size of the relativistic corrections at this nearly Newtonian compactness. The functions with I > 4 are smaller than those 
shown (in both panels) by a factor of order 10~^ or smaller and are not shown. 
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FIG. 6: The same mode as in Fig. ^, but for a strongly relativistic uniform density star (n — 0) with compactness M/R — 0.15. 
Upper left panel: The axial-parity fluid functions Ui{r) for / < 6. The Newtonian functions (dashed curves) are also shown for 
comparison. Upper right panel: The polar-parity fluid functions Wi{r) and Vi{r). Lower left panel: The axial metric functions 
hi{r) for I < 6. Lower right panel: The polar metric functions Hi^i{r). In all cases, the functions with / > 6 are of order 0.1% 
or smaller and therefore not shown. 
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FIG. 7: The same as Fig. H but for a relativistic n = 1 polytrope. 
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FIG. 8; A polar-led hybrid mode of a relativistic uniform density barotrope. The Newtonian counterpart to this mode (dashed 
curves) is the m = 2 polar-led hybrid with frequency kn ~ 1.232 (see Paper I). As in Fig. ^, the Newtonian and relativistic fluid 
functions are indistinguishable because the mode is shown in the weakly relativistic regime, with compactness M/R = 10"*. 
The left panel shows the fluid functions Ui{r), Vi{r) and Wi{r) while the right panel shows the metric functions hi{r) and 
Hi^i{r), all for / < 3. The functions with Z > 3 are smaller than those shown (in both panels) by a factor of order 10~* or 
smaller and are not shown. 
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FIG. 9: The same mode as in Fig. but for a strongly relativistic uniform density star (n = 0) with compactness M/R — 0.15. 
Upper left frame: The axial fluid functions Ui{r). The Newtonian functions (dashed curves) are also shown for comparison. 
Upper right frame: The polar fluid functions Wi{r) and Vi{r). Lower left frame: The axial metric functions hi(r). Lower right 
frame: The polar metric functions Hi^i{r). In all cases, the functions with I > 5 are of order 0.5% or smaller and therefore not 
shown. 
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FIG. 11: Gravitational radiation reaction timescales for the fastest growing I = m Newtonian r-modes (dashed hnes) and their 
relativistic axial- hybrid counterparts (solid curves). The timescales are shown as a function of compactness for uniform density 
stars (left panel) and n = 1 polytropes (right panel) of fixed baryon mass, Mb = IAMq. 
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FIG. 12: Gravitational radiation reaction timescales for the I = m = 2 Newtonian r-mode (dashed lines) and its relativistic 
counterpart (solid curves). The timescales are plotted versus compactness for fixed polytropic index (left panel) and versus 
polytropic index for fixed compactness (right panel). All of the stellar models have the same baryon mass, Mb = IAMq. 



